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■ Abstract 

o ■ 

^ | ' We present an analysis of QCD sum rules for pion form factor in next-to-leading order of 

perturbation theory for the case of pseudoscalar pion currents. The essential instanton con- 
tribution is reanalysed with account for present more accurate values of parameters entering 
Single Instanton Approximation (SIA). The theoretical curve obtained for Q 2 dependence of 
pion form factor is in a good agreement with existing experimental data. To calculate NLO 
^sD . corrections for double spectral densities we developed an effective computational technic. The 

details of the method together with the results for pion form factor in a more theoretically clean 
case of axial interpolating currents will be presented elsewhere. 
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! 1 Introduction 
Ok 

One of the first applications, following the discovery of the fact, that perturbation theory can 
be applied to study exclusive high momentum transfer processes HI E] , was the estimates 
of asymptotic Q 2 — > oo behavior of electromagnetic form factors of hadrons. It turned out, 
that for moderate momentum transfers the agreement between asymptotic QCD formulas and 
experimental data for pion and proton form factors is very poor. The latter is not surprising, if 
one recalls, that the validity of pQCD predictions depends in a large extent on the possibility 
to neglect power corrections 0(Q~ N ), where Q is the momentum transferred, and leave only 
leading contribution. For pion the leading l/Q 2 -contribution in the region Q 2 — > oo is given 
by hard rescattering of valence quarks (one gluon exchange subprocess). But, a straightforward 
use of the asymptotic formalism in the region Q 2 < 10 GeV 2 leads to the conclusion, that the 
mean virtuality of the exchanged gluon is smaller then 300 MeV. Such contributions should 
be accounted for as pion wave function corrections, which receive additional power suppression 
and at large Q 2 scale only as 1/Q 4 . Also, in this region finite size or nonperturbative effects 
become important and as a consequence, a complete evaluation of all contributions within 
pQCD approach utilizing light-cone wave functions at moderate values of momentum transfer 
is impossible. 

First systematic study of finite size corrections to asymptotic pQCD formulas was performed 
in [1J[3]. It was concluded, that an account for transverse momenta of quarks as well as Sudakov 
corrections makes pQCD predictions self-consistent at momentum transfers as low as few GeV. 
Moreover, it was argued that pQCD prediction dominated form factors even at moderate values 
of momentum transfer. However, a more careful analysis 0, specifically taking into account an 
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"intrinsic" pion transverse wave function, lead to the conclusion, that even self-consistency of 
pQCD predictions was achieved, it is still insufficient to describe experimental data. 

On the other hand, an alternative approach I10L II H I12j based on direct evaluation of 
form-factors from QCD sum rules easily accounts for experimentally measured values of pion 
electromagnetic form factor F 7T {Q 2 ). In this approach a leading order result at moderate values 
of Q 2 is given by soft wave function overlap mechanism and scales as l/Q 4 at large momentum 
transfers. Since then, the importance of inclusion both hard (pQCD) and soft contributions to 
describe experimental data was further emphasized and discussed. Later, in this paper, we will 
focus on QCD sum rule approach, which significantly reduces model dependence of the results 
obtained. So, let us discuss different contributions to pion electromagnetic form factor in the 
language of QCD sum rules. 

The method of QCD sum rules [Jj is designed to estimate low-energy characteristics of 
hadrons, such as masses, decay constants and form factors. Within this framework we can easily 
go into the region of intermediate momentum transfers and take due care about nonperturbative 
effects. 

It is already become a standard 1 , when considering the region of moderate momentum 
transfers for pion form factor, to account both for hard scattering and soft wave function overlap 
mechanism. Within QCD sum rule approach the soft contribution is dual to the lowest-order 
diagram, while the hard contribution is given by diagrams having higher order in coupling 
constant a s and as a consequence suppressed relative to soft contribution with additional factor 
a s /TT ~ 0.1. This extra suppression is in a complete agreement with asymptotic behavior of the 
pion electromagnetic form factor, calculated in pQCD 

F ha rd(Q2) = M^) f 1 dx f 1 dy ^M = ^1 (1) 

9 J J xyQ z Q 2 

At asymptotically high Q 2 the 0(a s /ii) suppression of hard contribution is more than com- 
pensated by its slower decrease with Q 2 . However, such compensation does not occur in the 
region of moderate momentum transfer, where, as we have seen, the soft contribution may be 
as important as hard contribution. 

The pion electromagnetic form factor was studied using a bunch of different frameworks, like 
QCD sum rules P EB HU H2] , light-cone sum rules [EJ HH US] , and different pQCD based ap- 
proaches employing convolution with phenomenological pion distribution amplitudes. It should 
be mentioned that within all these approaches correlation functions with axial interpolating cur- 
rents were studied. Another possibility is to consider a correlation functions with pseudoscalar 
currents, used as pion interpolating current,s|16| I17j. This simple change of pion interpolating 
current has however crucial consequences for the effects responsible for deviation of pion elec- 
tromagnetic form factor from pQCD predictions. It is known, that in scalar and pseudoscalar 
channels there are large instanton corrections due to direct interaction of quarks with instantons 
[TBI Ell- 
in this paper we consider NLO QCD sum rules for pion form factor, where pseudoscalar cur- 
rents are used as pion interpolating currents. As we already mentioned, in this case in contrast 
to the case of axial interpolating currents, there is large instanton contribution given by zero- 
modes of quark propagators in external (anti-)instanton field. Recently, a lot of quantitative 
self-consistent information about instanton properties in QCD was gained from phenomeno- 
logical estimates and numerical simulations. For a review of recent advances see |20j . As a 
consequences, today we are able to take relevant instanton contributions well under control. 
Another difference from the conventional approach employing axial currents is the asymptotic 
of pion electromagnetic form factor at large momentum transfers. In the pseudoscalar case 

1 See, for example, |S] 
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in limit Q 2 — > oo pion electromagnetic form factor scales as 1/Q 4 |2H 122 j . There is simple 
explanation of this fact 2 : the well known 1/Q 2 asymptotic corresponds to the scattering of 
virtual photon on the pion in twist 2 state, while taking pseudoscalar currents we artificially 
put pion in twist 3 state. From this discussion, we expect that perturbative contributions to 
pion electromagnetic form factor in pseudoscalar case should be of the same order as higher 
twist corrections in the axial case. Our analysis shows, that it is really the case. The obtained 
results of pion form factor calculation are in good agreement with existing experimental data 
as well as with predictions of sum rules based both on pseudoscalar and axial currents. There 
is however discrepancy in the values of different contributions compared to previous sum rule 
analysis based on pseudoscalar currents jTH] . The source of discrepancy is discussed in the main 
body of the paper. While being consistent with the different sum rule predictions, our results 
differ from the predictions made in |17j . 

Besides rich physical content of pseudoscalar channel, which draw our attention to this prob- 
lem, we had in mind the development of effective technic for calculation of QCD corrections to 
double spectral densities. So, it is natural, that we chose as a playground the simplest situation. 
In our forthcoming publications we are planning to present the details of our method as well as 
perform a detail analysis of NLO QCD sum rules with axial interpolating currents. Here, we 
would like to stress that an inclusion of radiative corrections is conceptually very important, as 
only in such a way we can simultaneously account for both hard and soft contributions. However 
in present case they turn out to be numerically small. 

The paper is organized as follows. In section 2 we describe our framework and give explicit 
expressions for next-to-leading order corrections to double spectral density. Next, in section 3 
we discuss instanton contributions, supplementing OPE of QCD sum rules in this case. Section 
4 contains our numerical analysis. Finally, in section 5 we draw our conclusions. 



2 Derivation of QCD sum rules 

To determine pion electromagnetic form factor we employ the approach of three-point QCD sum 
rules. This procedure is similar to that of two-point sum rules and the information from the 
latter on the pion coupling to its current is required in order to make predictions for the pion 
electromagnetic form factor. Here we are using pseudoscalar interpolating currents to describe 
pion states. The vacuum to pion transition matrix element of pseudoscalar current is defined 
by 

m 2 2 

where f n = 131 MeV and (0|^|0) | j Q e y = — (0.23GeV) 2 . Next, the pion electromagnetic 
form factor is given by hadronic matrix element of electromagnetic current: 

WWMp)) = U(Q 2 )(p, + p») + f-(Q 2 )(p, -Pp, (3) 

where j^ 1 = euUj^u + edd'jfj.d, the momenta of initial and final state pions were denoted by p, 
p' and Q 2 = —q 2 (q = p — p') is square of momentum transfer. Conservation of electromagnetic 
current leads to relation between two form factors introduced above. It is easy to obtain, that 

f-(Q 2 ) = U(Q 2 ) s - 1 Q^, (4) 

where the following notation was introduced: si = p 2 , S2 = p' 2 ■ 



2 We are grateful to A.Khodjamirian clarifying this point to us 
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Following the standard procedure for the evaluation of form factors in the framework of 
QCD sum rules, we consider the three-point correlation function: 

n M (p,j/,g) = i 2 J dxdye^'- x - ^\Q\T{u{x)^d{x),jf{<d),u{v)^d{y)) (5) 

Similar, to pion form factor decomposition U^(p,p' , q) is presented as a sum of two independent 
Lorentz structures: 

n M = n + (p M +^) + n_(p M - P ;). (6) 

The scalar amplitudes IT are the functions of kinematical invariants, i.e. IT = Ui(p 2 ,p' 2 , q 2 ). For 
the calculation of QCD expression for three-point correlator one employs the operator product 
expansion (OPE) for the T-ordered product of currents. As a result of OPE one obtains besides 
leading perturbative contribution also power corrections, given by vacuum QCD condensates. It 
turns out, that condensate contributions are small in our case, so in what follows we will neglect 
them. As, we already mentioned, in the case of pseudoscalar interpolating currents there are 
essential non-perturbative corrections due to instantons. In this section we will consider only 
perturbative contribution to sum rules, while the instanton contribution will be considered in 
next section. 




Figure V. LO diagram 



Next, for the perturbative contribution for each of two scalar amplitudes we write a double 
dispersion relation in variables s\ = p 2 and s 2 = p' 2 at q 2 < 0: 

nf er V,lA Q 2 ) = 7^T2 [ /^T 2 ' 9 ! ^ d32 + ^tractions (7) 
(2vr) 2 J (si - p 2 ){s 2 - p' z ) 

The integration region in (JJJ) is determined by condition 3 

-i< <i (8) 

X L ' 2 {si,s 2 ,q 2 ) 

and 

X(x 1 ,x 2 ,x 3 ) = (x 1 +x 2 - x 3 ) 2 - Ax x x 2 . (9) 

At leading order in coupling we have only one diagram depicted in Fig. 1, contributing to 
three-point correlation function. At next to leading order we have 6 diagrams shown in Fig. 2. 
The calculation of corresponding double spectral density was performed with the standard use 
of Cutkosky rules. In the kinematic region q 2 < 0, we are interested in, there is no problem in 



3 In the present case this inequality is satisfied identically. 
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Figure 2: NLO diagrams 



applying Cutkosky rules for determination of p%(si, s 2 , Q 2 ) and integration limits in s% and s 2 . 
The non-Landau type singularities, not accounted for by Cutkosky prescription, do not show 
up here. It is easy to find, that at Born level p+(si, s 2 , Q 2 ) is given by: 



py(s 1 ,s 2 ,Q ) = fc3/2 , (10) 

where k = A(si, s 2 , — Q 2 )- Calculation of radiative corrections to double spectral density is 
in principle straightforward. One just needs to consider all possible double cuts of diagrams, 
shown in Fig. 2. However, the presence of collinear and soft infrared divergences together with 
ultraviolet divergences related to renormalization of pseudoscalar currents calls for appropriate 
regularization of arising divergences at intermediate steps of calculation and makes the whole 
analytical calculation quite involved. In this paper we present only the final expression for 
radiative corrections to double spectral density, leaving all details of the calculation to separate 
publication. We checked, that all infrared divergences cancel, while the remained ultraviolet 
divergences are subtracted with renormalization constants of pseudoscalar currents. 

P+( S1 ,S 2 , Q 2 ) = ( Sl ,S 2 , Q 2 ) + (g) P { l\s U S2, Q 2 ), (11) 

LO spectral density was given before, and for NLO correction we have the following expres- 
sion: 



PY(si,S2,Q ) = 3p72"-{ 



-6sis 2 log 2 f -— ) - 12sis 2 log f -— ) + 12sis 2 log f -— ) log f - — 
+12 Sl S2 log ( -Hi) log (-^) - 12 Sl s 2 log ( 1 - -) log (- — 

V V2) V x 2j V z 2j V x 2 

-12 SlS2 log 2 (— ) - 12 SlS2 log 2 (-Hi) -6si S2 log 2 

\SlJ V s 2j V V2 

2 2 ( x 2 

+2ir sis 2 + 33sis 2 + 9s\Zi + 9s 2 zi + 9Q z\ — 18sis 2 log — 

V • s i 

+18sis 2 log ( -— ) + 12sis 2 log ( — ) log ( -— ) - 6sis 2 log ( - — 

V s 2j \SlJ V s 2j V V2, 
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-12*1*2 log ( ^ ) log ( ) + 12 S1S2 log ( ) log ( 1 - Z ± 
.SlJ \ yi) V 2/2/ V ^2 

l.5>- L .s 2 log ( |i-J - 18sis 2 log J + 18sis 2 log f-^J - 12sis 2 Li 2 ^ 

•12sis 2 Li 2 I --^ ) + 12sis 2 Li 2 ( — ] — 12sis 2 Li 2 ( -— ) - 12sis 2 Li 2 ( 



+12sis 2 Li 2 - 12sis 2 Li 2 f - 12sis 2 Li 2 - 12sis 2 Li 2 fl - J 

}, (12) 
where the following notation was introduced: 

si = ^i-s 2 -Q 2 )-^, (13) 
X2 = \(si-s 2 -Q 2 ) + ^Vk, (14) 



y2 2 ,-. . -v ■ 2 



yi = \(si + Q 2 -s 2 )-\yfk, (15) 

I( Sl + Q 2 - S2 ) + I^, (16) 

*i = ^(5i + s 2 + Q 2 )-^, (17) 

z 2 = ^(si + s 2 + Q 2 ) + ^ (18) 



In the limit Q — > oo our NLO double spectral density takes the following asymptotic form: 

= i^{-4„|L_ 40 -_ 28 + 261og (-)_ 281og (^} (19) 

So, in our case of pseudoscalar currents in contrast to the case of sum rules with axial interpolat- 
ing currents and light cone sum rules double logarithms cancel. This result may be considered 
strange at first. It is known, that Sudakov factor, suppressing transverse momenta is required to 
make pQCD predictions self-consistent at large momentum transfers 4 . However, as was already 
noted in Introduction, in pseudoscalar case taking limit Q 2 — » oo we are computing not the 
leading asymptotic, but higher twist or power correction to the leading asymptotic from axial 
case and thus there is no contradiction to what we would ordinary expect. Moreover, employing 
the same technic for the three-point correlation function with axial currents 5 , we obtain already 
known expression for double logarithmic terms. 

Now, let us proceed with the physical part of three-point sum rules. The connection to 
hadrons in the framework of QCD sum rules is obtained by matching the resulting QCD ex- 
pressions of current correlators with spectral representation, derived from a double dispersion 
relation at q 2 < 0: 

1 f p phys (si s 2 Q 2 ) 

Ui(pl,pl,q 2 ) = ——^ / — l - o- 1 ' 2 ' ^dsids 2 + subtractions. (20) 

(2ny J (si -p{){s 2 -pi) 

Assuming that the dispersion relation (12()f) is well convergent, the physical spectral functions 
are generally saturated by the lowest lying hadronic states plus a continuum starting at some 



4 See, for example (23| 

5 The axial case will be considered in a separate publication. 
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thresholds s\ and s 2 . 

pf hys ( Sl ,s 2 ,Q 2 ) = P r s (s 1 ,s 2 ,Q 2 ) + 

9( Sl - sf) ■ 8(s 2 - sf) ■ p™ n \s u s 2j Q 2 ), (21) 

where 

P r s (s 1 ,s 2 ,Q 2 ) = (0|d 7 5«k-(p , ))(^(P , )l/i(Q 2 )l^(p))(^(p)|S75d|0) • 

(2tt) 2 5(si)5(s 2 ) + higher state contributions (22) 

In our approximation of massless quarks we put m 2 = 0. The continuum of higher states is 
modeled by the perturbative absorptive part of IT, i.e. by pi. Then, the expressions for the form 
factors fi can be derived by equating the representations for three-point functions IT from l|7|) 
and 1)20(1 . This last step constitutes a formulation of QCD sum rules for our particular problem. 



3 Instanton contribution 

Now, let us supplement the conventional operator product expansion in QCD sum rules, pre- 
sented in the previous section, with direct instanton contributions. This contribution could 
be taken into account quite accurately using semiclassical approximation, in which all gauge 
configurations are replaced by an ensemble of topologically nontrivial fields: instantons and 
anti-instantons. Since in the sum rules the correlation function is being probed mainly at dis- 
tances x, y ~ 0.2 fm, which are smaller compared to instanton spacing R ~ 1 fm, it is natural to 
use for calculations single instanton approximation (SIA). In this approximation the correlation 
function is dominated by a single instanton, the closest one. The effects of other instantons are 
taken into account at the mean-field level. Another attractive feature of this approximation is 
the possibility to carry out all calculations analytically. It turns out, that the bulk properties 
of instanton ensemble could be very well described with simple parametrization |24j : 

n(p) = n5(p-p), (23) 

with the average (anti-)instanton density and size are given by 

1 1 

n = - fm~ 4 , p = - fm. (24) 

2 F 3 V ; 

The actual calculation of the closest instanton contribution could be performed with the help 
of quark propagator in an instanton background field: 

Sjix, y; z) = SP(x, y; z) + ST m {x, y; z). (25) 

The zero- mode part of the propagator is given by |23 [22l [22j 

1 ~li 



Si (x,y,z)- — 
where 



T T 



(x - z)4>(y - z), (26) 



and m stands for current quark mass. The corresponding expression in the anti-instanton field 
is obtained through substitutions 

1-75 1 + 75 - + fnQ , 
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In the limit of small distances, which is the case in the present analysis, the expression for 
nonzero-mode part of propagator is simplified and is given by 

Sr m (x,y,z)^S (x,y), (29) 

where Sq denotes the free quark propagator. Next, in SIA the effects of other instantons different 
from the closest one are taken care of by substituting the current quark mass m in the expression 
for zero-mode part of quark propagator with effective mass m*. In the simplest approximation, 
it can be extracted from the value of quark condensate |28| : 

2 

m* = m — -7r 2 p 2 (uu) . (30) 

Now, it is quite straightforward to calculate an instanton contribution to the three-point 
correlation function ©, we are interested in. It is given by 

ri™\p,p',q) = -^^T^y J dAx f dA y f dAz j*-**-** x 

y 2 (y - z)n - (y + z) 2 y^ 



+ 



(x 2 + p 2 ) 3 \y\(y 2 + p 2 ) 3/2 (y + z) 4 \z\(z 2 + p 2 fl 2 V^m ^ 



,(31) 



where the sum over instanton and anti-instanton contributions together with integration over 
their positions were performed. In what follows we will perform the numerical analysis for 
pion electromagnetic form factor using Borel scheme. A corresponding Borel transformed Borel 
instanton contribution was calculated in [16j . The resulting expression could be found in the 
section with numerical results after the definition of Borel transformation. 



4 Numerical analysis 

For numerical analysis we are using so called Borel scheme of QCD sum rules. That is, to get 
rid of unknown subtraction terms in Q let us exploit the Borel transformation procedure in 
two variables s\ and Si- We define the Borel transform of three-point function IIj(si, S2, q 2 ) as 



MM 2 ,M 2 ,q 2 



B 12 Tli(si,S2,q 2 



lim 

n,m— >oo 



nl 



n „m+l 
H 

ml 



_d_ 

dsi 



sx=mMl ,S2=nM% _ 



Ui(si,s 2 ,q 2 



Borel transformation (|32j) of Q and (|2U|) gives 

1 



<S>i{MlMi,q') 



(2vr) 2 



oo roo 

ds\ I ds2 exp 
Jo 



Si 

M 2 



S2 

M 2 



(32) 



(33) 



In the following we put M 2 = M| = M 2 . If M 2 is chosen to be of order 1 GeV 2 , then the 
right hand side of ()33|) in the case of physical spectral density will be dominated by the lowest 
hadronic state contribution, while the higher state contribution will be suppressed. 

Equating Borel transformed theoretical and physical parts of QCD sum rules we get 



L(0)(M 2 ,Q 2 ) + ^1 



a s (M)J 



a s (M) 1142 



4tt 

+ ^a> ( MW)-4^ 

47r m* z (p) 



81 



•0(°)(M 2 ,Q 2 



(f ms \p 2 Q 2 ,p 2 M z ) 



4 

I 2 



^(OI^IO) 2 -^^ 2 ) 



(34) 
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where F n (Q 2 ) = f + (Q 2 ), p n is the normalization point [p, n ~ 0.5 GeV). We also introduced the 
following notation: 



(o) (M 2 ,Q 2 ) 
(1) (M 2 ,Q 2 ) 



(2vr) 2 
1 

Wf 



dx exp 



dx exp 



x 

w 

X 

w 



dyp+\s 1 ,s 2 ,Q 2 ), 
dyp ( ± ) (s 1 ,s 2 ,Q 2 ), 



(35) 



where x = s\ + s 2 and y = s\ — s 2 . Here, for continuum subtraction we used so called "tri- 
angle" model. To verify the stability of our results with respect to choice of continuum model 
we checked, that the usual "square" model gives similar predictions for pion electromagnetic 
form factor provided so ~ 1.5si is chosen. For more information about different continuum 
subtraction models see [TU] , 
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0.1- 
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Figure 3: Q 2 dependence of pion electromagnetic form factor 



The instanton contribution c/> mst (p 2 Q 2 , p 2 M 2 ) has the following expression d6: 

i 

aj3 



r s \zi,z 2 ) 



da 



d(3e~ azi e"^ +a) e *U-/s*a) • 



H{uo)H{a) 
2af3 



a + (3 



A 



z 2 

(2a(3 -A)- 



az\ + 
a 3 (3 



Az 2 A 2 



A±{1 - (3z 2 ) 
(3z 2 z 2 - 8(1 - 0z 2 ) r 
16 (l-(3z 2 ) 2 
aP(a + (5) 



z 2 A 



(36) 



where A = + a/5, uj = a = and H(z) = I\(z) — Iq{z) is defined in terms of the 
modified Bessel functions I n (z). 

In the above sum rules for pion electromagnetic form factor we also have taken into account 
the evolution of pion pseudoscalar current J(p) from scale M down to normalization scale p n : 

J(p 2 ) = J{pi)U(pi,p 2 ), (37) 

where 



+ 



a s (pi) - a s (p 2 ) 7i 
4?r 0i 



~72 


/V 


) 


.71 


ft. 





(38) 



Q 2 F n (Q 2 ) 
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Figure 4: Q 2 dependence of perturbative contribution to pion electromagnetic form factor 



and 



ih 
71 

72 



2 11- -N 



19 

- 1 ! 51 - -N f 



-2(-3C F ), 



(39) 

(40) 
(41) 
(42) 



Here C A = N c , C F 



2N C ' x 'c 



N c - number of colors and Nf denotes number of active flavors. 



FJ1 GeV 2 




Figure 5: Borel mass M 2 dependence of pion electromagnetic form factor at Q 2 = 1 GeV 2 

For the continuum threshold we use the value sq = 3.6 GeV 2 extracted from corresponding 
two-point sum rules for pion [22]. The relatively large separation of continuum from lowest 
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resonance in pseudoscalar channel (compared to vector and axial channels) is explained by 
almost complete dominance of the pion at low energies in this particular channel. A few words 
should be said about the value of m*(p) used in numerical analysis. The analysis of ^1 uses 
the value extracted from the quark condensate Ea. (|30j) . Later, a systematic analysis of SIA 
30 showed, that in the instanton vacuum with realistic density the square of effective mass 
parameter m* 2 (p) is not given by the square of effective mass m*(p) extracted from quark 
condensate Eq. (|30j) . The real value of m* 2 (p) is sufficiently smaller. In our analysis we use 
m* 2 (p) = (70 MeV) 2 . Next, we use two- loop renormalization group running of strong coupling 
constant with Aq CD = 325 MeV. In Fig. [21 we plotted the dependence of full electromagnetic 
pion form factor at fixed Q 2 = 1 GeV 2 (sum of instanton and perturbative contributions) 
on the value of Borel parameter M 2 . We see that "stability plateau" starts to develop for 
M 2 > 1.2 GeV 2 . To determine Q 2 -dependence of pion form factor in what follows we fix 
M 2 = 1.2 GeV 2 . As a result we get the results for pion electromagnetic form factor shown 
in FigEl (solid line is the sum of instanton and perturbative contributions, curve with long 
dashes denotes instanton contribution and curve with short dashes stands for perturbative 
contributions.). For completeness, in Fig^Jwe plotted the relative magnitude of LO and NLO 
perturbative contributions solid line is the sum of LO and NLO, curve with long dashes denotes 
LO contribution and curve with short dashes stands for NLO contributions.). 

We see, that obtained results for pion electromagnetic form factor are in good agreement with 
available experimental data. Our results are also in agreement with the predictions of sum rules 
employing axial pion currents. The dominant contributions to pion electromagnetic form factor 
is given by instantons. The difference in relative sizes of perturbative and instanton corrections 
computed here compared to similar quantities from |16j is explained by wrong value of m* 2 (p) 
used in the latter analysis. On the other hand, the small value of perturbative contribution is in 
agreement with our expectations: it should correspond to higher twist correction to leading twist 
2 contribution extracted from sum rules with axial interpolating currents for pion. Also, our 
perturbative correction has the same value of magnitude as similar higher twist contributions 
considered in the framework of light-cone sum rules in |14] 115]. So, we may conclude here that 
pion electromagnetic form factor is almost saturated by instanton contribution in pseudoscalar 
channel. It should be mentioned, that uncertainties of theoretical predictions in pseudoscalar 
case are still large (conservatively it is up to 50%) and are mainly due to strong dependence of 
results on the values of normalization point p n , Aq CD and square of effective mass m* 2 (p). These 
uncertainties could be further reduced by considering simultaneously sum rules for two-point 
correlation function of pseudoscalar currents. 

5 Conclusion 

We presented the results for pion electromagnetic form factor in the framework of three-point 
NLO QCD sum rules with pseudoscalar interpolating currents for pions. We redone the analysis 
of instanton contribution to three-point correlation function in this case. It was found, that while 
overall prediction for pion electromagnetic form factor made in ^B] holds true, the relative 
magnitude of different contributions is completely different. The instanton contribution almost 
saturates the prediction for pion form factor. The theoretical curve obtained for Q 2 dependence 
of pion form factor is in a good agreement with existing experimental data. Here, we for the first 
time computed radiative corrections to three-point sum rules with pseudoscalar currents. While, 
these corrections turned out to be small numerically the technic developed and tested here will 
certainly have numerous applications in determining QCD corrections to the form factors of 
other light mesons. We are planning to present in our forthcoming publications both the details 
of NLO calculation together with the results for pion form factor in a more theoretically clean 
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case of axial interpolating currents. 
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